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Abstract 

We provide numerical constructions of one-dimensional hyperuniform many-particle distribu- 
tions that exhibit unusual clustering and asymptotic local number density fluctuations growing 
more slowly than the volume of an observation window but faster than the surface area. Hyperuni- 
formity, defined by vanishing infinite-wavelength local density fluctuations, provides a quantitative 
metric of global order within a many-particle configuration and signals the onset of an "inverted" 
critical point in which the direct correlation function becomes long-ranged. By targeting a specified 
form of the structure factor at small wavenumbers using collective density variables, we are able to 
tailor the form of asymptotic local density fluctuations while simultaneously measuring the effect of 
imposing weak and strong constraints on the available degrees of freedom within the system. This 
procedure is equivalent to finding the (possibly disordered) classical ground state of an interacting 
many-particle system with up to four-body interactions. Even in one dimension, the long-range ef- 
fective interactions induce clustering and nontrivial phase transitions in the resulting ground-state 
configurations. We provide an analytical connection between the fraction of contrained degrees 
of freedom within the system and the disorder-order phase transition for a class of target struc- 
ture factors by examining the realizability of the constrained contribution to the pair correlation 
function. Our results explicitly demonstrate that disordered hyperuniform many-particle ground 
states, and therefore also point distributions, with substantial clustering can be constructed. We 
directly relate the local coordination structure of our point patterns to the distribution of the void 
space external to the particles, and we provide a scaling argument for the configurational entropy 
of the systems when only a small fraction of the degrees of freedom are constrained. By empha- 
sizing the intimate connection between geometrical constraints on the particle distribution and 
structural regularity, our work has direct implications for higher-dimensional systems, including 
an understanding of the appearance of hyper uniformity and quasi-long-range pair correlations in 
maximally random strictly jammed packings of hard spheres. 
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I. INTRODUCTION 



The relationship between the local structure of a many-particle system and interparticle 
correlations is fundamental to condensed-matter theory. This intimate connection provides 
a useful image of the regularity jl] of all phases of matter, allowing researchers to track 
the local structure over increasing length scales approaching the global system. In practice, 
one measures pair correlations between distinct points in the form of the structure factor 
S(k), which is proportional to the scattering intensity from x-ray or small-angle neutron 
scattering ^j. It is intuitive from such measurements that a hierarchy of structural order 
can be established, ranging from crystalline structures such as Bravais lattices js| to highly 
disordered systems, the prototypical example of which is the ideal gas ^-7|. Unfortunately, 
quantitative descriptors consistent with this stratification of order are difficult to identify, 
and this area of research is currently open. One recently introduced order metric 8] involves 
the notion of hyperuniformity of point patterns, whereby infinite-wavelength local density 
fluctuations vanish (4, 6]. This order metric explicitly indicates the degree to which density 
fluctuations are suppressed on large length scales. 

The local structure of a hyperuniform many-particle configuration (i.e., on the order of 
a few nearest-neighbor distances between particles) is by definition indicative of the global 
arrangement of particles 4|. Also known as superhomogeneity J9], this phenomenon is 
fundamental to the description of all Bravais lattices, lattices with a multiparticle basis, 
quasicrystals, and certain disordered systems possessing pair correlation functions decaying 
to unity exponentially fast [6] . We emphasize that while hyperuniformity in periodic config- 
urations is a trivial consequence of their intrinsic long-range order, the fact that disordered 
many-particle systems can also display this property is nonintuitive. This behavior is espe- 
cially surprising since the appearance of hyperuniformity marks the onset of an "inverted" 
critical point in which the structure factor vanishes in the limit of small wavenumbers while 
the direct correlation function, defined through the Ornstein-Zernike formalism, becomes 
long-ranged [l]. 

Hyperuniform systems have played a fundamental role in our understanding and design 
of materials, including those with large, complete photonic band gaps [lcj, "stealth" ma- 
terials invisible to certain frequencies of radiation and prototypical glassy structures 



consisting of maximally random strictly jammed (MR J) monodisperse hard spheres 
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Other ex amp les of disordered hyperuniform systems include noninteracting spin-polarized 



fermions |14| . Il5| , the ground state of liquid helium 16| , the density fluctuations of the early 
Universe |17{|, one-component plasmas and so-called g 2 - mvar iant processes ji], in which 
the form of the pair correlation function is held fixed over a certain density interval. Note for 
equilibrium many-particle configurations at positive temperature, hyper uniformity implies 
that the isothermal compressibility vanishes; this relationship does not hold, however, for 
nonequilibrium systems. 

Hyperuniform particle distributions possess structure factors with a small-wavenumber 
scaling S(k) ~ k a for a > 0, including the special case a = +00 for periodic crystals. 
This behavior implies that the variance u 2 N {R) in the number of particles within a local 
observation window (here a d- dimensional sphere of radius R) increases asymptotically as 
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a%(R) 



R d ~ l In R, 

R d ~ a , 

R d -\ 



a = 1 
a < 1 
a > 1 



(R +00). 



(1) 



however, all known hyperuniform configurations to date have a scaling parameter a > 1 
HQ , meaning that the second asymptotic regime of the number variance in ([1]) has never 
been observed in either theoretical or experimental studies. Indeed, the aforementioned 
MR J packings, which are maximally disordered among all jammed sphere packings with 
diverging elastic moduli, possess a small-wavenumber scaling a = 1, and this observation 
has provoked the question of whether this value corresponds to a minimal scaling among all 
hyperuniform point patterns. Zachary, Jiao, and Torquato have provided strong arguments 
that this claim is indeed true for strictly jammed hard-particle packings [ljj, but it is unclear 
whether general point patterns must also possess exponents a > 1. Here we provide for 
the first time constructions of "anomalous" disordered hyperuniform many-particle ground 
states for which a < 1, demonstrating the diversity of possible structures within of this class 
of systems. 

Our approach involves placing explicit constraints on the so-called collective coordinates 
associated with a point distribution, which are defined b y a Fourier transform of the local 
density variable (discussed in Section II below) 



21] . Controls on collective coor 



18 



dinates have been previously used in the development of novel stealth materials 



ll| and 



in the identification of unusual disordered classical ground states for certain classes of pair 



potentials 22|. This problem can be viewed as the determination of the ground state of a 
many-particle system with up to four-body interactions [18||; duality relations that relate 
the energy per particle of a many-body potential in real space to the corresponding energy 
of the dual (Fourier-transform ed) potential can be used to examine analytically the ground 
state structures and energies |23j. Importantly, since collective coordinates directly probe 
the configuration space associated with the two-particle information of the structure fac- 
tor, they are ideally suited to the construction of hyperuniform point patterns. Formally, 
we numerically construct a configuration of particles whose spatial distribution is consis- 
tent with a targeted form of the structure factor at small wavenumbers. By constraining 
a certain number of degrees of freedom in the system, we "fix" the positions of a known 
fraction of the total number of particles based on the locations of the remaining particles 
and the implicit constraints imposed by the targeted form of S(k). By varying the fraction 
of constrained degrees of freedom within the system, we are able to explore directly the 
relationship between hyper uniformity and internal structural constraints of a many-particle 
configuration, allowing us to interpolate between the "disordered" and "ordered" regimes of 
hyperuniformity. 

In order to elucidate the connection between the local coordination structure and pair 
correlations for our anomalous hyperuniform ground states, we have investigated the distri- 
bution of the available void space external to the particles. Prior work on MRJ packings 
of binary hard disks has shown that the appearance of hyperuniformity in a many-partic 



system is related to the underlying distribution of the local voids between particles [13]; 
in this sense, the void space is more fundamental to the local structure than the particles 
themselves. Strong arguments have also been put forth to support the claim that exponen- 
tial values a less than unity in the small-wavenumber region of the structure factor indicate 
the presence of larger interparticle voids with higher 



requency, thereby deregularizing the 
microstructure while maintaining hyperuniformity [13j|. This behavior is notable since it is 
not obvious that hyperuniformity can be consistent with a highly clustered microstructure; 
see Fig. [TJ Here we provide further evidence to link rigorously the void space and the local 
coordination structure of a point pattern, and we highlight the differences in the void space 
distribution for "regular" and "anomalous" hyperuniform systems. Since we can directly 
control the fraction of constrained degrees of freedom via collective coordinates, our results 



FIG. 1. (Color online) Numerically- generated configurations of particles in two dimensions with a 
circular local observation window of radius R. Both configurations exhibit strong local clustering 
of points and possess a highly irregular local structure; however, the configuration on the left is 
hyperuniform while the one on the right is not. The hyperuniform point pattern was generated 
with the same methodology outlined in Section III of the text. 

have implications for understanding how the void space distribution is affected by increased 
constraints on the many-particle configuration. Indeed, our work directly supports the fun- 
damental role of the void space in the microstructure and reinforces the relationship between 
constraints on the local structure and the aforementioned observed minimal scaling a = 1 
found in S(k) for MRJ hard sphere packings. 
Our major results are summarized as follows: 

(i) Disordered hyperuniform many-particle ground states can, counterintuitively, exhibit 
a substantial degree of clustering in the absence of a large number of constraints on 
the particle distribution (Sections IV and V). 

(ii) The order-disorder phase transition that occurs upon increasing the fraction of con- 
strained degrees of freedom is related to the realizability of the constrained contribution 
to the pair correlation function (denned below) (Section IV). 

(iii) Hyperuniform particle distributions with anomalous asymptotic local density fluctua- 
tions (i.e., slower than the volume but faster than the surface area of an observation 



window) can be constructed, and these fluctuations are intimately related to the dis- 
tribution of void sizes external to the particles (Section V). 

(iv) With few constrained degrees of freedom (e.g., a perturbation from an ideal gas), the 
entropy (configurational degeneracy) decreases linearly with the number of constraints 
imposed on the particle distribution (Section V). 

Section II provides a brief overview of the important ideas related to point processes, 
collective coordinates, and hyperuniformity. We apply these concepts in Section III to discuss 
how control over collective coordinates can be used to numerically generate configurations 
of hyperuniform point patterns, including those with anomalous asymptotic local density 
fluctuations, to a high numerical precision. Section IV explores how increasing the fraction 
of constrained degrees of freedom within hyperuniform systems affects the observed pair 
correlations and, therefore, the local coordination structure. In Section V we provide explicit 
calculations for the void statistics of our hyperuniform point patterns under weak and strong 
constraints, and we draw explicit connections among the regularity of the local structure, the 
exponential form of the small-wavenumber region of the structure factor, and the distribution 
of the local voids. Concluding remarks are given in Section VI. 

II. STOCHASTIC POINT PATTERNS, COLLECTIVE COORDINATES, AND 
HYPERUNIFORMITY 

We consider many-particle configurations to be realizations of stochastic point processes 
in some subset of Euclidean space M, d . A (finite) stochastic point pattern is formally defined 
as a distribution of N points {r^} in some compact space V of volume (Lebesgue measure) 
V. We consider the case where the distribution is statistically homogeneous with periodic 
boundary conditions on V; the thermodynamic limit N, V — > +00 with p = N/V = constant 
can be taken appropriately to extend the point pattern to Euclidean space M. d . The statistics 
of the process are determined by an iV-particle probability density function P^(r N ), which 
need not be a Gibbs measure. Equivalently, one can specify the countable set of generic 
n-particle probability density functions p n (r n ), defined by 
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The function p n is therefore the probability density associated with finding a subset of 
any n particles within volume elements dr n . Note that for statistically homogeneous point 
patterns p\ = p. Related to the generic n-particle probability density function is the n- 
particle correlation function g n (r n ), defined by 

P n g n {v n )= Pn {v n ). (3) 

Of particular importance is the pair correlation function <72( r ), which can be made integrable 
by subtracting its long-range value of unity to give the total correlation function h(r) = 
g 2 (r) — 1. A Fourier representation of g 2 (r) is given by the structure factor S(k), defined by 

S(k) = l + ph(k), (4) 

where we utilize the following convention for the Fourier transform: 

/>)= / exp(-ik.r)/(r)dr. (5) 

Corresponding to any single configuration of points {r N } is a local density variable 

JV 

p(r) = ^(r-r;), (6) 

where 5 denotes the Dirac delta function. The ensemble average of this local density with 
respect to the statistics of the point process is 

(p(r)} = P, (7) 

and the autocorrelation function is given by 

(p(r 1 )p(r 2 ))=pS(r)+p 2 g 2 (r) (8) 

with r = r 1 —r 2 . Note from ([8]) that the autocorrelation function contains two contributions: 
a delta function corresponding to the self-correlation of a point in the process and the 
pair correlation function between two distinct particles. The self-correlation contribution is 
independent of the distribution of particles in the system and arises for all correlated and 
uncorrelated point patterns. 

For statistically homogeneous point patterns subject to periodic boundary conditions, it 
is convenient to assume ergodicity and equate ensemble averages with volume averages over 
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the unit cell. This assumption is expected to be valid in the thermodynamic limit. One can 
show that the volume-averaged local density and autocorrelation function are 



p(r) = p (9) 
1 N 

p(x + r)p(x) = p5(T) + - 5 ( r - r ^)> (10) 

where = Vj — r^. Equation ffTUj) suggests the following alternative definition of the pair 
correlation function: 



3^=1 

Since the Dirac delta functions in ( TTU1) are by definition localized, this result has little 
practical utility when handling finite particle distributions. However, one can take advantage 
of the periodicity of the unit cell to expand the local density in a Fourier series according to 

1 N 
j=i k 

which is equivalent to a discrete (inverse) Fourier transform. The wavevectors k in (TT2l) are 
determined by the geometry of the unit cell; if the unit cell is formed with basis vectors 
{ej}, then the wavevectors satisfy 

k ■ ej = 2ixm (13) 

for all i and for some m G Z. For simplicity, we will henceforth consider a d- dimensional 
cubic cell [0,L] d C M. d , which implies k = 2iin/L for some n e Z d . Rewriting (pD in the 
form 

k 

where 

N 

P( k ) = 5^ ex P( _zk ' r i)' ( 15 ) 

3=1 

we observe that the local density is the discrete (inverse) Fourier transform of p, which we 
call a collective density variable. 

The identity ()9]) can also be obtained using the Fourier representation (|T2|) . meaning that 
only the mode k = contributes to the local density on average. However, the autocorrela- 



tion function is now of the form 

p(x + r)p(x) = {pY1 ex P( ik • r ) + Tn Yl Yl exp ^ k ' ( r ~ r ^ ^ 



k i+i k 

+ E ex p(^ k • r ) Op( k )i 2 - ^] . ( 17 ) 

k 



which, by comparing with ([8]), implies [24] 

92(r) = ^- 2 J2^^-r)[\p(k)\ 2 -N}. (18) 

k 

The result f|T8|) allows one to directly compute the pair correlation function from the collec- 
tive density variables p; note that the k = mode must be included in this calculation to 
ensure the correct long-range behavior g 2 {r) — > 1 as ||r|| — > +oo. In practice, one must trun- 
cate the wavevector summation in ( JT8l) . leading to oscillatory approximations to g 2 within 
some threshold determined by the cut-off magnitude of the wavevectors. 

Hyperuniform point patterns constitute a subclass of point processes lacking infinite- 
wavelength local density fluctuations 4]. Specifically, it has been shown that the variance 
<J%{R) in the number of points within a local spherical observation window W(R) of radius 
R and volume v (R) scales asymptotically as |4j 



a%(R) = (N(R)) [A N (R) + B N (R)/R + lower-order terms] , (19) 

where (N(R)) = pv(R) is the average number of points in the observation window. The 
coefficients A^(R) and B^(R) in ( !T9|) are determined solely by the two-particle information 
of the point pattern: 



A N (R) = l + p h(r)dr (R -> +oo) (20) 

Jw(R) 

B "W = - wU/$m L M > ryr * ( ^ +oo) - (21) 



So long as h(r) — > faster than r d , the leading-order coefficient An(R) converges asymp 



totically as A^(R) = A^ = lim|| k ||^ 5(k) [25]. By definition, a hyperuniform point pat- 



tern possesses a number variance growing slower than the volume v(R) of the observation 
window (equivalently, the mean number of points (N(R))), implying that Ajq = and 
infinite-wavelength density fluctuations vanish. 

The most common examples, including all Bravais lattices, periodic non-Bravais lat- 
tices, quasicrystals possessing Bragg peaks, and certain disordered point patterns with pair 

10 



correlation functions decaying to unity exponentially fast, of hyperuniform point patterns 
possess constant number variance coefficients Bn(R) = [4]. This behavior implies that 
the isotropic structure factor S(k) possesses a small-wavenumber scaling Dk a with a > 2, 
including the special case a = +00 for periodic structures. However, it is also possible to 
find hyperuniform point patterns for which < a < 2, in which case C\ < B^(R) < C2R 
as R — > +00 for some constants C% and C2. The most well-known examples of these types 
of "anomalous" local density fluctuations occur when S(k) ~ k as k — > 0, in which case 
Bn{R) = A\ \n(R) + A2 with A\ and A2 constant. This situation has been well-characterized 
in three-dimensional maximally random jammed packings of hard spheres 12], the ground 
states of liquid helium [if]], and noninteracting spin-polarized fermion ground states jl^j ]. 
However, examples where a < 1 have heretofore not appeared in the literature. 



III. COLLECTIVE COORDINATE CONSTRUCTION OF HYPERUNIFORM 
POINT PATTERNS 

One goal of this work is to construct examples of hyperuniform point patterns possessing 
the aforementioned "anomalous" asymptotic local density fluctuations, meaning that the 
number variance grows slower than the volume of an observation window but faster than 
the surface area. Collective density variables provide an attractive means to control the 
small-wavenumber region of the structure factor S(k), thereby allowing us to construct a 
hyperuniform point pattern with targeted local density fluctuations. Specifically, we define 
an objective function $ according to 

Hr N ) = J2[S^,r N )-S (k)} 2 , (22) 
keQ 

where S (k) is the targeted form of the structure factor and Q denotes some finite subset 
of wavevectors k. The structure factor is determined using collective density variables; 
specifically, 

S(k;r N )= l -^l (MO), (23) 

where p(k), implicitly a function of the particle positions r^, is defined by ( I15p . The zero- 
wavevector is excluded from (|23|) since it provides an O(N) contribution to the structure 
factor, corresponding to a delta function in the thermodynamic limit from the long-range 
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behavior of g>2- By expanding (|22|) . one can show that our minimization problem corre- 
sponds to finding the classical ground state of a many-particle system with up to four-body 
interactions 

§(r N ) = ^ Vi(Ti,Tj,Tt,T m )+ ^ V 3 (T i ,Tj,Tt)+^r i V2(Ti,T j )+Vo, (24) 

where 

v a {y u Tj, r e , r m ) = — ^ cos ( k • r ij) cos ( k ' Y tm) (25) 

keQ 

4 x - 

v 3 (r h Tj, r i> = jpZ^ cos ( k • r ij) cos ( k • r u) (26) 
keQ 

^( r - r i) = |^cos(k- ri ,)[l-5o(k)] (27) 
keQ 

^ = ^[5 (k)-l] 2 . (28) 

keQ 

The set Q in ( 12 2 p is chosen to contain all wavevectors, excluding the zero mode, with 
norm less than some upper bound K. This construction allows us to target specifically the 
small- wavenumber region of the structure factor, which controls the asymptotic local density 
fluctuations. The target function So is chosen with the form 

5 (k) = D||k|| a forallkGQ. (29) 

In order for the target function to correspond to a realizable point pattern, it is necessary 
that D > to enforce positivity of the structure factor. The parameter a determines the 
asymptotic behaviors of the pair correlation function and the number variance [c.f. (|IJ)]. 
Previous work [18| has considered the cases a = 1,2,4,6,8, and 10 in dimensions d = 2 
and 3. It has recently been conjectured that a = 1 corresponds to the minimal exponent 
consistent with the constraints of saturation and strict jamming in sphere packings 



131; 



however, systems for which a < 1 have not been reported in the literature, and their 
statistical properties are unknown. The objective function f )22|) is minimized to within 10~ 17 



of its global minimum using the MINOP algorithm 26|, |27j, which has several c omp utational 
advantages for this type of investigation as previously reported in the literature 18j ■ MINOP 
applies a dogleg strategy that uses a gradient direction when one is far from the minimum, 
a quasi-Newton direction when one is close, and a linear combination of the two when one 
is at intermediate distances from the minimum. 
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It is important for this study to verify that the constructed point patterns are indeed 
hyperuniform with the correct targeted asymptotic local density fluctuations. This criterion 
requires high resolution of the small-wavenumber region of the structure factor. Specifically, 
the smallest observable wavenumber magnitude in the collective coordinates representation 
(in a d- dimensional cubic unit cell) is k mm = 2n/L = 2np 1 / d /N 1 / d , where L is the box length, 
N is the number of particles, and p is the number density. To ensure hyperuniformity, we 
therefore require that lim7v->+oo S(k m in) = 0, where the limit is taken at constant density. 

Since any simulation necessarily requires choosing N finite, it is essential to select a value 
of N sufficiently large to enforce both hyperuniformity and the desired form of the structure 
factor near the origin. Unfortunately, the 0(N~ 1 / d ) scaling of /c m ; n makes obtaining such 
resolution increasingly difficult in higher dimensions. Our interest is in verifying the existence 
of anomalous hyperuniform point patterns and understanding their statistical properties, 
and we therefore limit our studies to one dimension, where the scaling is most favorable, with 
N = 2000 particles. It should be appreciated, however, that hyperuniform point patterns 
with logarithmically-growing asymptotic density fluctuations are known in arbitrarily high 
dimensions [14J. Importantly, since our minimization procedure is equivalent to finding the 
classical ground state of a long-range interaction with up to four-body potentials and can 
be used in principle to construct hyperuniform point patterns in any dimension, nontrivial 
phase behaviors can still be observed 28j , and we are therefore able to extend our conclusions 
to higher-dimensional structures. 



IV. COLLECTIVE COORDINATES AND REALIZ ABILITY OF POINT PAT- 
TERNS 

For a general d- dimensional point pattern of N particles, there are dN translational de- 
grees of freedom in the absence of constraints on the system. One must therefore choose a set 
of wavevectors Q for the objective function ( 1221) containing only a fraction x of these degrees 
of freedom. In one dimension there are 2M(K) = floor (KL/n) wavevectors, excluding the 
zero mode, with magnitude less than or equal to K. Inversion-invariance of the modulus of 
the collective density variable implies that M(K) of these wavevectors can be independently 
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constrained; we therefore define a new parameter 

which represents the fraction of independently constrained degrees of freedom from the 
objective function $. 

For the case where the targeted structure factor <So(k) = for all k e Q, it has been 
previously shown 2^] that increasing the parameter x induces a greater degree of order on 
the particle distribution. Specifically, in one dimension the corresponding point patterns 
are disordered for < x < 1/3 an d crystalline for x > 1/2 {^J; intermediate values of % 
interpolate between these two regimes [3^]. However, it is known that target functions of 
the form (|2"§|) interfere with this order- disorder phase transition; here we provide analytic 
results suggesting that this transition is shifted to higher values of \ f° r all finite a 31]. 

For a one-dimensional point pattern, the wavevectors are of the form k = 2irm/L for 
m G Z, and one can write the collective density variable as: 

N 

p(m) = exp(— Vlitmrjl L) . (31) 

Additionally, the total correlation function is of the form [cf. ( [18]) ] 

2 + °° 



h ^ = E cos ( 27rmr / L ) [\p(^m/L)\ 2 - N] , (32) 

m=l 

which for the targeted point pattern can be decomposed as: 



M 2 +°° 

h(r) = — ^cos(2?rmr/L) [D(2mn/L) a - 1] +— ^ cos(2vrmr/L) [|p(2vrm/L)| 2 - iV] 

m=l m=M+l 

(33) 

= /i (r;M) + / il (r;M), (34) 

where ho(r;M) is the contribution to the total correlation function due to constrained 
wavevectors and hi(r; M) is the unconstrained contribution. The function ho can be simpli- 
fied as 

h (r; M) = ( - + N l a D ) cos(2vrmr/L)m a - A ^ cos(2vrmr/L) (35) 

^ ' m=l m=l 

M 

= C(a, D) cos(27rmr / L)m a - (2/iV) cos[(M + l)nr/L] csc(7rr/L) sin(M7rr/L) 

m=l 

(36) 
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where 

2 a+l ir a D 

C (a,D) = —- (37) 

is a parameter-dependent constant. The global minimum of ho(r; M) occurs at r = 0, 
corresponding to 

M 

h (0;M) = C(a,D)J2 m<x ~ ( 2M / N ) ( 38 ) 

m=l 

= C(a, D)H { - a \M) - (2M/N), (39) 

where 

n 

H (a \n) = J2 m ~ a ( 40 ) 

m=l 

is the harmonic number of order a. 

The negative contribution to h (0; M) in (13"9~|) suggests that there may be an upper 
threshold M* beyond which ho(0; M) < 0. For any values of M in this region, the constrained 
contribution h to the total correlation function of the point pattern is no longer in itsel 



realizable as a point process. The realizability problem in classical statistical mechanics 



22| 



and the associated iV-representability problem in quantum statistics [33J are notoriously 
difficult and unsolved problems in physics that ask under what sufficient and necessary 
conditions a reduced two-particle correlation function can be expressed as the integral over 
a full iV-particle probability density. In the classical case, one can consider specifying a pair 
correlation function g 2 and attempting to construct a corresponding point process. Known 
necessary realizability conditions on g 2 include 

g 2 (r) > for all r (41) 
S{k) > for all k (42) 

along with the somewhat weaker Yamada condition 

a 2 N (R)>9(l-9) (43) 



on the fractional part 9 of the average number of particles in an observation window 34] 



The Yamada condition appears easy to satisfy in all but relatively low dimensions 32| . The 



determination of other realizability conditions on g 2 is an open problem 351 ] . 

Figures [2] and [3] compare the pair correlation functions and the constrained contributions 
h (r) + 1 for numerically-constructed point patterns (using the methodology of Section 
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FIG. 2. (Color online) Left panel: Pair correlation function g 2 for numerically-constructed hype- 
runiform point patterns with small-wave number scalings Dk a and \ = 0.1. Right panel: Con- 
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FIG. 3. (Color online) Left panel: Pair correlation function g 2 for numerically-constructed hype- 
runiform point patterns with small-wavenumber scalings Dk a and x = 0.35. Right panel: Con- 
strained contributions to the pair correlation functions. 



Ill) with small-wavenumber exponents a = 0.5,1.0, and 2.0 and x = 0-1 an d 0.35. For 
X = 0.1, corresponding to a small fraction of constrained degrees of freedom, the constrained 
contribution h (r) + 1 places only moderate constraints on the local structure of the system, 
primarily controlling oscillations in g 2 beyond approximately five nearest-neighbor distances. 
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FIG. 4. Typical threshold values x*( a ) beyond which the constrained contribution ho(r) to the 
total correlation function is no longer realizable as a point process. This curve corresponds to 
choosing S(K) = 0.5, where K is the magnitude of the maximally constrained wavevector. Note 



that as a — > +oo we recover the crystallization threshold x* = 0-5 reported in Ref. 



2C|. 



Interestingly, the small-r behaviors of g2{r) and h (r) + 1 are strikingly different. Although 
the constrained contribution to the pair correlation function generates an effective repulsion 
between particle pairs, the full pair correlation function indicates a tendency for particles 
to cluster at short pair separations. It follows that the unconstrained contribution to the 
pair correlation function plays a substantial role in determining the local structure for this 
system. 

However, the situation is quite different upon increasing the constrained degrees of free- 
dom to x — 0.35. Figure [3] shows that the constrained contribution to #2 almost exactly 
mirrors the full pair correlation function, implying that sufficiently constraining the collec- 
tive density variables places a strong constraint on the local structure of the point pattern. 
It follows that the value M* beyond which h Q (0; M) < is an indicative precursor to the loss 
of realizability of the targeted structure factor. We have mapped the threshold value M* 
(equivalently, x*) 111 Fig. HI We emphasize that this loss of realizability is associated with 
negativity of the real-space pair correlation function; the structure factor itself is still positive 
over its entire domain. Interestingly, as the exponent a controlling the small-wavenumber 
region of the structure factor increases, we recover the value x = 0-5 corresponding to crys- 
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tallization in the case where 5o(k) = for all k G Q. This observation suggests that the 
threshold values of \ beyond which /io(0; M) < generalize this phase transition. In Section 
V, we provide additional arguments to support this claim. 



V. VOID STATISTICS AND COORDINATION STRUCTURE 
A. Exclusion probability functions 

The n-particle correlation functions contain information concerning the relative locations 
of points within a point process, and, in principle, specifying the countably infinite set (in 
the thermodynamic limit) of such functions is sufficient to completely determine the point 
pattern. However, any finite collection of correlation functions contains only partial details 
of the spatial arrangements of the points, implying that there are degenerate structures 



with these same statistics 



361 ] . In particular, the n-particle correlations functions do not in 



themselves provide direct information about the space exterior to the points, or the so-called 
void space. It has been shown for pent patter Q, H (and r a„do m .nedia ft that the 



distribution of the void space is indeed a more fundamental descriptor of the point process 
than the arrangements of the points themselves. Here we are interested in characterizing 
the relationship between asymptotic local number density fluctuations and the void space 
statistics; in particular, we would like to examine the constraints that the exponent a in 
the small-wavenumber region of the structure factor places on the distribution of the void 
space. 

One can define two types of "exclusion" functions, both of which measure the availability 
of empty space surrounding points of a stochastic process. The void exclusion probability 
function Ey(r) is the probability of finding a d- dimensional spherical cavity of radius r 
centered at an arbitrary position in M d . The void exclusion probability has recently been 

role in the covering and quantizer problems from discrete 



shown to play a fundamental 
geometry and number theory 39[ . Closely related to this descriptor is the particle exclusion 
probability function Ep(r), which is the probability of finding a d- dimensional sphere of 
radius r centered on a point of the point process but containing no other points. Figure [5] 
highlights the differences between these functions. 

The exclusion probability functions are complementary cumulative distributions of the 
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FIG. 5. Events contributing to the void exclusion probability Ey{r) (upper left) and the particle 
exclusion probability Ep(r) (lower). The points correspond to a realization of a disordered point 
process. 



void and particle nearest-neighbor functions Hy(r) and Hp(r), respectively 37|, |38(. The 
void nearest-neighbor function is the probability density of a finding the nearest point of a 
point process with respect to an arbitrary location in IR d within a radial distance r + dr. The 
particle nearest-neighbor function is defined similarly but with respect to nearest neighbors 
between two points of a point process. One therefore has the following simple relationships 
between these sets of functions: 



H v (r) 
H P (r) 



dEvir) 

Or 
dE P (r) 

dr 



(44) 
(45) 



One can relate the void and particle exclusion probabilities via a simple probabilistic 
construction [15j. Specifically, we consider a generalized exclusion probability Ey(r;e), 



which is the probability of finding a d- dimensional annulus of outer radius r and inner 
radius e; by definition, Ey(r; 0) = Ey(r). Taking the derivative of this function with respect 
to the inner radius e gives a function proportional to the probability of finding a point within 
a small radial region inside the annulus and the annulus itself devoid of points. It follows 



19 



that Ep(r), the conditional probability of finding a spherical cavity centered on a point, is 



£ P r = lim — , 46 

e^o+ ps(e) oe 

where s(e) is the surface area of a <i-dimensional sphere of radius e. This construction is 
known in the theory of point processes [40!] and has also been used in the literature to identify 
the void statistics of certain point patterns related to problems in number theory, random 
matrix theory, and quantum mechanics [14J]. One can without loss of generality define the 
exclusion correlation function 77 (r) according to 

"W - Wr (47) 

This function provides a measure of the correlations between neighboring points in a stochas- 
tic point pattern and is identically unity for a Poisson point process. It is interesting to note 
that for a system of equilibrium hard spheres of diameter D, the exclusion correlation func- 



tion is given by [3 

\[E v {r)Y\ r<D 
v{r) = C ~ (48) 

[[E v (D)}-\ r>D, 

which depends only on knowledge of Ey(r) and is monotonically nondecreasing for all r with 
77(0) = 1. 

Further insight into the probabilistic meanings of Ep, Ey, and r\ can be gained by intro- 
ducing the notion of the particle space, defined to be the subset (of Lesbesgue measure zero) 
of M. d occupied by the points of the point process. The particle exclusion probability function 
Ep(r) is then the fraction of the particle space that can be decorated by a d- dimensional 
sphere of radius r containing no other points of the process. To define the void exclusion 
probability Ey(r), one decorates all of the points in the process by spheres of radius r and 
then determines the fraction of all space not occupied by the spheres; this value corresponds 
to the portion of space available to insert a cavity of radius r [39( . The exclusion correlation 
function 77 (r) then provides a measure of the relative available space for a cavity of radius r 
in the particle space compared to the external void space. 

Torquato and coworkers 



371 ] have provided the following series representations for the 
20 



exclusion probability functions: 



Mr) = l + ^Y(kTT) /^(r fe )n^(l|x-r,||;r)dr, (49) 
k=l ^ ' J 3=1 

+00 ( \k p ^+1 

Mr) = 1 + E r(k + T) J 9k ^ rk+1) II m(l|ri " 'iM** ( 5 °) 

k=\ 3= 2 

where m(r; i?) = Q(R—r). Since these functions are special cases of a more general canonical 
n-particle correlation function [41| . one can establish rigorous upper and lower bounds by 
truncating these series at finite order. Specifically, by writing 

+00 

Ev/p(r) = T, E v/p(r), (51) 



fc=0 



where Ey) p = 1, we have the following hierarchy of bounds: 



E v /p(r) < £ E$ p (r) (£ even) (52) 

k=0 

i 

E v/P (r)>Y^E^ P (r) (I odd), (53) 



fc=0 

which become sharper with increasing I. 



B. Local statistics of anomalous hyperuniform point patterns 

We have been able to successfully construct point configurations exhibiting anomalous 
asymptotic local number density fluctuations. Figure [6] provides images of the structure 
factors for our configurations at % — 0.1 and \ = 0.35. As we show, for all wavevectors 
within the constrained portion of the spectrum, the structure factor matches its target value 
within an exceedingly small numerical tolerance (on the order of 10 -17 ). In addition to 
the systems shown, we have also been able to reliably construct configurations with small- 
k exponential behaviors a > 0.25. In order to keep the exposition clear, we have only 
presented results for a = 0.5 with the disclaimer that our conclusions will apply for other 
point patterns with anomalous local number density fluctuations. 

It is interesting to note the substantial differences in the structure factors of the systems 
for x — 0.1 and x = 0.35, particularly for unconstrained wavevectors. For x = 0.1, the 
structure factor exhibits an unusually slow decay to its asymptotic value of unity; we have 
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FIG. 6. (Color online) Left panel: Structure factor with small-A; behavior Dk a (inset) for 
numerically-constructed hyperunifbrm point patterns with \ = 0-1- Right panel: Structure factor 
with small-A; behavior for x = 0.35. 

fit the large- A; region of the structure factor with an asymptotic fit of the form 1 + ft/k 1 
and have found a power-law decay 7 = 1. This behavior is due to the local clustering of 
particles as expected from the small-r region of the pair correlation function in Fig. [2j This 
effect can be directly observed in Figs. [7] and which provide illustrative portions of our 
numerically-constructed hyperuniform point patterns at x — 0-1 an d X = 0.35. As has been 



previously reported in the literature 



11 



18], increasing the fraction of constrained degrees 
of freedom in the many-particle system has the effect of imposing greater local order in the 
form of an effective short-range repulsive interaction. By increasing x from 0.1 to 0.35, the 
relative influence of the constrained wavevectors on the pair correlation function increases, 
suppressing the formation of local clusters. However, we also observe that as the exponent 
a controlling the small-wavenumber region of the structure factor decreases (equivalently, 
as anomalous local number density fluctuations appear), this effective repulsion between 
particles becomes noticeable weaker, manifested in the pair correlation function by larger 
values of #2(0). This behavior suggests that anomalous hyperuniform point patterns possess 
greater variability in their local structures, particularly with regard to the shapes and sizes 
of voids between particles. 

We have verified the expected asymptotic behaviors of the number variance for our 
numerically-constructed point patterns as shown in Fig. |9j The asymptotic scalings of 
these fluctuations exactly correspond to their theoretical predictions. In particular, we 

22 



FIG. 7. (Color online) Portions of numerically-constructed hyperuniform point patterns with 
X = 0.1 and small-/c exponential scalings (upper left) a = 0.5, (upper right) a = 1.0, and (lower) 
a = 2.0. 
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FIG. 8. (Color online) Portions of numerically-constructed hyperuniform point patterns with 
X = 0.35 and small- A; exponential scalings (upper left) a = 0.5, (upper right) a = 1.0, and 
(lower) a = 2.0. 

have provided the first example of a hyperuniform point pattern for which the asymptotic 
number variance grows more slowly than the volume of an observation window but faster 
than a logarithmic scaling. Interestingly, the local clustering of points at x = 0.1 gen- 
erates strong oscillations in the number variance that persist for several nearest-neighbor 
distances. In contrast, these local oscillations essentially vanish after two nearest-neighbor 
distances at x — 0.35, reflecting the strong constraints placed on the local structure by the 
small-wavenumber region of the structure factor. 

Our calculations for the void and particle exclusion probabilities of these systems, shown 
in Figure [TOj demonstrate previously-unobserved statistical properties for hyperuniform 
point patterns. For a Poisson point pattern, one has the result that 

E v (r) = E P (r) = exp[-pv(r)], (54) 

implying that the exclusion correlation function r)(r) = 1 for all r. This result follows from 
the absence of interparticle correlations for the process and the underlying Poisson probabil- 
ity distribution for the number of particles within an arbitrary compact set. Gabrielli and 
Torquato [42| have provided strong arguments to suggest that for any hyperuniform point 
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FIG. 9. (Color online) Left panel: Number variances afj(R) for numerically-constructed hyperuni- 
form point patterns with % = 0.1. Right panel: Corresponding number variances for \ = 0.35. 

pattern, the void exclusion probability should asymptotically decay faster than for a Poisson 
point process. This behavior implies that arbitrarily large cavities within the system, while 
not prohibited by the constraint of hyperuniformity, are expected to be significantly rare 
events owing to the underlying regularity of the global structure of the pattern. It is there- 
fore not unreasonable to expect the functional form of Ey(r) for the Poisson point process to 
provide an upper bound on the exlcusion probability of any hyperuniform point pattern, and 
this observation is indeed rigorously true for point patterns generated from fermionic parti- 
cle distributions (so-called determinantal point processes) 141 ] . More generally, the Poisson 
result will place an upper bound on Ey for any point pattern with n-particle correlation 
functions g n < 1 for all n. 

For x — 0.1, we observe the unusual property that Ey(r) is greater than the Poisson 
result for all values of r that can be reliably determined from numerical simulation. It is 
instead the particle exclusion probability function that is bounded from above by the Poisson 
curve. To understand this discrepancy, we first note that Ey(r) and Ep(r) are rigorously 
bounded from below by [c.f. ( |52i) ] fl^ | 



Ey(r) > 1 — pv(r) 
E P {r) > 1 - Z{r), 



(55) 
(56) 



where Z(r) is the cumulative coordination number 

Z(R) = p I e(R -r)g 2 {r)dr. 



(57) 
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FIG. 10. (Color online) Left panel: Void- and particle-exclusion probabilities for the numerically- 
constructed hyperuniform point patterns (x = 0.1) along with the reference curve for the Poisson 
point process. Right panel: Corresponding functions for x = 0.35. 

These bounds become sharp at low density or small r. Therefore, while Ey(r) is related 
to the geometry of a cavity within the void space, the particle exclusion probability Ep 
depends explicitly on the local coordination structure of the underlying point process. 

To elucidate further the relationship between the local coordination structure and the 
void statistics, we can consider a modification of the number variance problem, whereby one 
measures fluctuations in the number of points within an observation window centered on a 
point of the point process. Let N^p (R) denote this quantity; it can be represented as 

N^(R)= J^Q(R- \\Ti-rS, (58) 
i=i 

where the prime on the summation means that particle i is excluded. The average value of 
this random variable is 

(N®(R)) = p [ g 2 (r)e(R - \\r\\)dr = Z(R). (59) 
Jm. d 

The cumulative coordination number therefore measures the local number density within 
the particle space. It follows that when Z(R) < pv(R), (N$(R)), the average number of 
points in the particle space within an observation window of radius R, is less than or equal 
to (N V (R)), the average number of points of the process within a window in the void space. 
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This behavior then implies that the points are more greatly dispersed within the particle 
space, and Ep(R) > Ey(R)\ equivalently, r)(R) > 1. Note that this analysis is consistent 
with the lower bounds ( 15 5 p and (I56p on the exclusion probability functions. 

Conversely, for the case where Z(R) > pv(R), we have that (nP(R)) > (N V (R)), which 
suggests that the points are more closely located within the particle space, leaving larger 
cavities within the void space. We therefore conclude that Ep(R) < Ey(R) [r](R) < 1], 
and the point process should exhibit local clustering among points. These claims are also 
consistent with our results for the exclusion probability functions in Fig. [10l whereby we 
observe a transition from 1](R) < 1 at x = 0.1 to r)(r) > 1 for \ = 0.35. Indeed, the 
configurations at x — 0-1 exhibit substantial clustering among points (c.f. Fig. [7]). 

Our arguments can be extended by examining the scaling of the configurational degener- 
acy with the fraction of constrained degrees of freedom \. Here we measure this degeneracy 
by calculating the entropy (logarithm of the degeneracy) of the system relative to an ideal 
gas of iV particles in a volume V on the line. For the ideal gas, we coarse-grain the system 
by dividing the volume V into M 3> 1 cells such that no more than one particle occupies 
each cell with probability one, thereby representing the degeneracy associated with the dN 
translational degrees of freedom as a combinatorial occupancy problem. The size of a cell 
determines the length scale, meaning without loss of generality that we need only consider 
the regime p = N/M <C 1. Assuming that the particles are indistinguishable, the number 
of configurations Q available to the system is 

a = (M-N)W: (60) 

Since the underlying distribution of particles is uniform within the cells, Boltzmann's formula 
for the entropy gives (with kp = 1) 

Ml 



S = \nn = \n , , (61) 

[(M-N)\N\\ , y 1 

which for large M and iV becomes 

S = M\nM — M — (M — N) ln(M — N) + (M — N) — N\nN + N (62) 

= -M ln(l - N/M) +N\n(M/N) + N\n(l - N/M). (63) 

Under the assumption that N/M < 1, we have the following result for the entropy per 
particle S ideal = S ideal /N: 

Sidcai = 1 - N/M + \n(M/N) w 1 - ln(p). (64) 
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The entropy of the ideal gas is therefore large and positive as expected. For the density 
regime p <C 1 that we have in mind, one can simplify further by taking S ~ — ln(p), which 
diverges to +00 for small p. 

Suppose now that we constrain K degrees of freedom, where K ^ N. This construction 
correponds to making a small perturbation away from the ideal gas configuration with N — K 
degrees of freedom still available to the many-particle system. We again divide the volume 
V into M cells of unit length (N/M <C 1). For sufficiently small values of x (i-e., near the 
ideal gas), we may assume that the length scale of the effective repulsion between particles 
is negligible compared to the cell size, meaning that the N — K unconstrained particles may 
be distributed freely among the M cells. Note that the K constrained degrees of freedom 
are explicitly determined once these particles have been placed. The number of microstates 
available to the system is then 

Ml 

^ = (M — N + K)\(N — K)\ ' (65) 
and the configurational entropy is 

\M-N + KJ \M - N + K J \ N-K J v ; 

Defining the fraction of constrained degrees of freedom x — K/N, we may write for the 
entropy per particle S — S/N: 

S = (l- M/N - x) ln(l - N/M + X N/M) + (1 - X ) ln(M/7V) - (I - x ) ln(l - X ) (67) 

where X <?C 1. For N/M <^ 1, this expression simplifies as 

S = l- X -(1-X) ln(l - X) + (1 - X) HM/N) - (N/M)(l - X f ■ (68) 

The last term in this expression is negligible within the density regime where Sidcai ~ 
ln(M/iV) = ln(l/p). Since ^idcai is large and positive compared to the first terms of this 
result, we have the expected scaling 

^— ~ 1 - X, (69) 

>J ideal 

suggesting a roughly linear decrease in the entropy of the system for small values of x- 

By increasing the parameter a, we expect also to increase the rate at which S — > with 
respect to x since higher values of a are associated with larger effective radii around the 
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constrained particles. However, since we require that our configurations are hyperuniform, 
there are additional implicit constraints on the "unconstrained" degrees of freedom. For- 
mally, hyperuniformity requires that the local structure of a point process approach the 
global structure over sufficiently short length scales, on the order of several nearest-neighbor 
distances. This behavior is typically associated with a highly regularized distribution of 
points, such as with a Bravais lattice. However, our results for Ey at x — 0.1 suggest an 
alternative mechanism by which hyperuniformity can be achieved in a point pattern. Specif- 
ically, the fact that i](r) < 1 is consistent with local clusters of particles that are globally 
regularized by an increased probability of finding sufficiently large voids to separate them. 
In this case, the appearance of these large voids external to the particle space is apparently 
essential to enforce hyperuniformity of the point pattern by overcoming the highly inho- 
mogeneous local structure of the clusters. In the context of our analysis above, for small 
perturbations from the ideal gas, the high configurational degeneracy that remains after 
constraining only a few degrees of freedom implies that highly disordered configurations are 
most likely to appear from our numerical constructions. However, with the added implicit 
constraint of hyperuniformity the system will sacrifice local structural regularity for cluster- 
ing of points that are globally separated by sufficiently large voids, resulting in a negatively 
correlated exclusion correlation function. 

Figure (TT] highlights this behavior by examining the cumulative coordination numbers 
Z(R) at x — 0-1 an d X — 0.35 for our systems. Although at x — 0-1 the a = 2 system 
exhibits greater clustering at small-r than for a = 0.5, at large r this trend reverses, which 
is consistent with the appearance of larger voids with higher probability and supports our 
claim that these voids serve to regularize the global structures of the systems. Using the 
lower bound ( 1561) on the particle exclusion probability function, we also observe explicitly 
the effect of the locally clustered structure on the small-r particle exclusion probability. 

Upon reaching x — 0.35, we recover the usual behavior associated with hyperuniformity. 
By constraining a sufficient number of degrees of freedom, the effective interparticle repul- 
sions induced by the collective coordinate constraints control the small-r region of the pair 
correlation function (and, therefore, the cumulative coordination number), prohibiting local 
cluster formation. The regularizing factor in this case is therefore the distribution of voids 
within the particle space, contained in E P {r). Indeed, the void exclusion probability E v {r) 
is highly constrained by this effective repulsion and decays to zero faster than Ep(r), result- 
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FIG. 11. Cumulative coordination numbers Z{R) for numerically-constructed hyperuniform point 
patterns with small- wave number exponents a = 0.5 and a = 2.0. The fractions of constrained 
degrees of freedom are x = 0.1 and x = 0.35. 



ing in a positively correlated exclusion correlation function. By increasing the exponent a 
governing the small-wavenumber region of the structure factor, we observe increased regu- 
larity in the local structure, corresponding to a decreased Z(R) and Ey(r) and a particle 
exclusion correlation function Ep(r) that decays more rapidly to zero, in perfect accordance 



with the a: 
Torquato 



orementioned void-space criterion on hyperuniformity put forth by Gabrielli and 



421 ] . It is also noteworthy that increased void-space constraints associated with 



increased a are consistent with the behavior of the structure factor for MRJ hard-sphere 
packings, about which we have more to say in Section VI. 



VI. CONCLUDING REMARKS AND DISCUSSION 

We have provided the first known constructions of disordered hyperuniform many-particle 
ground states possessing anomalous local density fluctuations. Such systems are defined 
by a number variance afj(R) that asymptotically scales faster than the surface area of 
an observation window but slower than the window volume. By controlling the collective 
density variables associated with the underlying point pattern, we have also be able to probe 
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the relationship between interparticle correlations and constraints on the local coordination 
structure. Specifically, we have provided detailed statistics to measure the distribution of 
the void space external to the particles, including measurements of the void and particle 
exclusion probabilities. 

Under sufficiently low constraints on the system, our numerically constructed many- 
particle distributions exhibit substantial clustering, resulting in a highly inhomogeneous 
local structure. However, on the global scale of the system as measured by asymptotic local 
density fluctuations, these local clusters are separated by comparatively large interparticle 
voids, thereby regularizing the microstructure and preserving the constraint of hyperuni- 
formity that we impose. Indeed, this effect becomes more pronounced upon passing from 
the "anomalous" regime of hyperuniformity to the more usual case, where o~ 2 N (R) ~ R d ~ l 
asymptotically (i.e., for all periodic point patterns, quasicrystals with Bragg peaks, and dis- 
ordered systems with pair correlations decaying exponentially fast) j^, 6]. Upon increasing 
the fraction of constrained degrees of freedom within the system, we are able to preclude this 
clustering affect by reinforcing the effective interparticle repulsion imposed by our targeted 
structure factor S(k). Furthermore, we have shown that this effective repulsion becomes 
increasingly more pronounced as the exponent a governing the small- wavenumber scaling of 
S(k) is increased. 

It follows from these observations that one can formally define an effective repulsive ra- 
dius around each point within a hyperuniform point pattern. However, so long as this radius 
is not substantially large compared to the expected interparticle spacing p -1 ^, i.e., in the 
absence of microstructural constraints, clustering effects can still dominate the interparticle 
correlations and the local coordination structure. We have shown that this effect is en- 
tropically favorable since slight deviations from the ideal gas are still associated with an 
exponentially large configurational degeneracy. However, this degeneracy is expected to in- 
crease rapidly upon constraining a sufficient number of degrees of freedom, or, equivalently, 
increasing the effective repulsive radius. We have shown that this loss of configurational 
degeneracy is associated with a highly-constrained void space distribution, which can be 
considered as a signature of predominantly "repulsive" hyperuniform point patterns. 

Our results have particular implications for understanding the appearance and nature of 
hyperuniformity in MRJ packings of hard spheres. It is interesting to note that equilibrium 
distributions of hard spheres are known not to be hyperuniform (4] except at the close- 
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packed density, at which point the system freezes into a crystal with long-range order. MRJ 
hard-sphere packings are therefore unique in that they are nonequilibrium systems that are 
uniformly mechanically rigid, and it is this rigidity that has been shown to be essential 
'or the onset of hyperuniformity (with logarithmic asymptotic local density fluctuations) 
1^ | . Importantly, rigidity places severe geometric constraints on the local arrangements of 
particles 43J and has been shown to regularize the void space distribution on the global 
scale of the microstructure We have demonstrated in this work that by decreasing the 
exponential form of the structure factor within the small-wavenumber region, these void- 
space constraints are relaxed in accordance with the decreased effective radius surrounding 
the particles. Since MRJ packings are maximally disordered among all strictly jammed 
packings, it follows that such systems must already possess the maximal number of degrees 
of freedom consistent with the geometric constraints of strict jamming. Therefore, any 
increase in the distribution of the void sizes is inconsistent with these same constraints, 
highlighting why exponents a < 1 in the small-wavenumber scaling of S(k) have never been 
observed for such systems. 



Our work has also raised a number of interesting questions related to the physics of 
collective coordinate constraints. The mathematical properties associated with collective 
coordinates are surprisingly subtle and have only partially been explored in the literature 
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2l|. Constraining a collective density variable, including, for example, either a complete 
suppression to zero or fixing its magnitude, results in a highly nonlinear equation relating 
the components of the particle positions {r.,}. For higher- dimensional systems, it has been 



previously observed 22j that these nonlinear equations have a tendency to "couple" in 
such a way that one must go beyond x = 0-5 to crystallize the system, meaning that one 
cannot simply count constraints and degrees of freedom. It is an open problem to determine 
analytically the relationship between x an d the "true" number of constrained degrees of 
freedom; recent work has examined the fraction of normal modes with vanishing frequency 
as a more appropriate indicator of the latter j^]. We have provided some analysis here 
in one dimension to suggest how the configuration space is constrained with increasing x 
by determining the entropy for small deviations from the ideal gas. Certainly for higher 
values of x our simple linear scaling will break down; characterizing the deviations from this 
linear behavior, especially in higher dimensions, is an attractive problem warranting further 
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consideration. 
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